Strong decays of heavy baryons in Bethe-Salpeter 

formalism 



Xin-Heng Guo* Ke-Wei Wei^ and Xing-Hua Wu^ 

Institute of Low Energy Nuclear Physics, Beijing Normal University^ 

Beijing 100875, China 

October 8, 2007 



Abstract 

In this paper we study the properties of diquarks (composed of u and/or d 
quarks) in the Bethe-Salpeter formahsm under the covariant instantaneous ap- 
proximation. We calculate their BS wave functions and study their effective 
interaction with the pion. Using the effective coupling constant among the di- 
quarks and the pion, in the heavy quark limit mq ^ cxd, we calculate the decay 
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1 Introduction 



A light baryon composed of u, d, or s quarks is a very complex three-body system in 
which all the three light quarks play important roles in the dynamics of the baryon. 
However, for a baryon containing a heavy quark, in the heavy quark limit, dynamics in 
the baryon is greatly simplified theoretically. In the heavy quark limit, the heavy quark 
effective theory (HQET) [1] shows that the light quarks move in an effective static color 
field (in the rest frame of the heavy baryon) and can not see the spin and flavor degrees 
of the heavy quark. The heavy baryon system has an extra SU (2) f^SU (2)^ symmetry. 
Therefore, the spin and isospin of the light quarks and the heavy quark are conserved 
separately. 

In recent years, more and more experimental results about heavy baryons have 
been reported by various experimental collaborations, e.g. the discovery and measure- 
ments of [2,3], [4], and HjJ" [5], etc. in addition to much more results about 
baryons containing one c or 6 quark [2]. However, the structures and properties of 
these baryons are not very clear, hence more precise experimental measurements and 
detailed theoretical studies are urgent. 

For a baryon with one heavy quark, the diquark picture has been taken in various 
references, see e.g. Refs. [6,7]. Since the spin and the isospin of two light quarks 
within the heavy baryon are conserved, they can be regarded as a two-quark system, 
"diquark", and the diquark then combines with the heavy quark to form the heavy 
baryon. In this diquark picture, the heavy baryon can be reduced to a two-body 
system. In this paper, we will focus on the study of heavy baryons containing one 
heavy quark and two light quarks {u and/or d). We shall not consider the effects of 
the isospin symmetry violation. 

The purpose of this paper is to calculate the decay widths of the strong decay 

processes Tjq^ — > Aq + vr with Q = c or b . These processes can be shown schematically 

in Fig. [TJ From this figure, we can see that the light diquark from the mother baryon 

TiQ combines with the spectator heavy quark to form the daughter baryon Aq after 

emitting a very soft pion. 

From Fig. [H we can see that, to calculate the amplitude (AQ(/^)7r(g) [^[^^(Ps)), 

(*) 

we need to know the wave functions of Eg and Aq and the effective interaction among 
the diquarks and the pion. In the Bethe-Salpeter (BS) equation approach, one can 
show that the baryonic wave functions appearing in these decay amplitudes are the 
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Figure 1: The schematic picture for the decay processes, Sq"* — > Aq + tt. The hght 
diquark from Eg'' combines with the spectator heavy quark to form Aq after emitting 
a very soft pion. 

BS wave functions of Tiq^ and Aq in the heavy-quark and hght-diquark picture, which 
have been given in Refs. [6,7]. On the other hand, since the pion emitted by the 
diquark is very soft, the interaction among the diquarks and the pion can be treated 
by chiral perturbation theory. In order to calculate the (lowest order) effective coupling 
constant among the diquarks and the pion, G-,,^^, we will establish the BS equations 
for the scalar diquark and the axial-vector diquark, respectively. Then Gt^ip^ can be 
expressed as the overlap integral of the BS wave functions of the scalar diquark and the 
axial- vector diquark. To simplify the BS equations for the diquarks to tractable forms, 
we will impose the so-called covariant instantaneous approximation in the kernels of 
these BS equations. This approximation was also applied to the BS equations for Sq'* 
and Aq [6,7]. Furthermore, we will assume that the kernels contain a scalar confinement 
term and a one-gluon-exchange term. Throughout our calculations, we will take the 
heavy quark limit (that is to say, to neglect all the l/rriQ corrections except for those 
coming from the kinematic factors). 

Now let us give some more detailed explanation about the points mentioned above. 
First, the decay amplitude shown in Fig. [T]can be written out in terms of the following 
equation: 

(^Q(^A)7r(g)|s[j*^(PE)) = j d^ixiX2yiy2Uv)x^^{x2,xi)SQ{xi - yi)"^Xp^ ,a(2/i, 1/2) 
X A-/ix2 - u)A-''-^\v - y2) J2(^iq)\T <P\umvm , (1) 



where Pa, Ps, and q are the momenta of A, S>, , and vr, respectively, Xp and Xp a are 
the Bethe-Salpeter (BS) wave functions of Aq and Sg-*, respectively, A,^ and A^-^ are 
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the propagators of the scalar diquark and the axial-vector diquark yj, respectively, 
Sq is the propagator of the heavy quark Q, and i is the color index. 

Second, let us discuss the effective coupling among the diquarks and the pion. Since 
the pion emitted by the diquark is very soft, one can expand the low energy effective 
Lagrangian of the diquarks and the pion in terms of the pion's momentum. To lowest 
order, we only need to consider the following interaction vertex: 

b,i 

where b is the isospin index. Then, the matrix element (7r(g)|T0*(M)^^(f)|O) in Eq. 
([T]) can be calculated out from this interaction vertex. 

On the other hand, the effective coupling constant can be related to a transition am- 
plitude (more details will be shown in Sect. [3]), 

^.,^^-^^(0), (4) 

Jtt 

where A^q"^) ~ q^{P4„ 'i\A~{q)\Pip, r, i), which can be presented by the overlap integral of 
the BS wave functions of the diquarks (p and (p. To calculate this transition amplitude, 
we need the BS wave functions of and ip. 

The diquark systems have been studied in Ref. [8] in which the authors first solve 
the BS equations in the rest frame of diquarks and then boost their solutions to a 
moving frame. Different from the approach of Ref. [8] , we will solve the BS equations 
directly under the so-called covariant instantaneous approximation in a general coor- 
dinate system. Furthermore, the kernel in our paper is different from that used in 
Ref. [8]: there are only Coulomb part (arising from one-gluon-exchange) and the scalar 
linear part in our kernel, V ~ + nr. 

One can see from Eq. ([T]) that, to calculate the decay amplitudes of heavy baryons, 
we also need the BS wave functions of heavy baryons. The BS wave functions of Aq 

and Sg^ have been studied in Refs. [6,7]. We will take the results given there as the 

(*) 

input to calculate the decay amplitudes of Sg in this paper. 

We will constrain the ranges of parameters in our model, (or m^p, which is 
related to m,^) and k^, by comparing the theoretical and experimental results about 
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the average momentum of the heavy quark in heavy baryons [10]. With these ranges 
of parameters, we calculate the decay widths of Eg Aq + vr. 

The remainder of this paper is organized as follows. In Sect. |21 we discuss the BS 
formalism under the covariant instantaneous approximation in some details. After a 
brief discussion about the BS equation of a general two-quark system, we derive the 
BS equations for the scalar diquark (p and the axial-vector diquark respectively. We 
also give the normalization conditions of the BS wave functions in this section. In Sect. 
El we calculate the effective coupling constant among the diquarks and the pion. In 
Sect, m the calculation of the decay widths of Sg Aq + vr are carried out in the 
BS formalism under the covariant instantaneous approximation. In Sect. O we discuss 
how to constrain the parameters in our model and give numerical results for the decay 
widths of Tjq^ Aq + vr. Sect. IHlis reserved for our conclusions and some discussions. 
We also include three appendices (jAl [HI and O) in this paper. Appendix |X] contains 
some definitions. We review some results of Refs. [6,7] in Appendix [Bl In appendix ICl 
we discuss the normalization conditions of the BS wave functions of heavy baryons. 



2 Bethe-Salpeter equations for diquarks 

In this section, we will derive the BS equations for a scalar diquark (p and an axial- 
vector diquark (f. Before doing this, we will first discuss the general BS formalism for 
a two-quark system (other than the system composed of a quark and an anti-quark). 
The BS wave functions of the two-quark system are defined as the following: 

xM^x,)fp = (0|TC(a:i)^^(x2)|P,A;), (5) 
xM,x,)'^: = (P,A;|TV^^(x2)*C(xi)*|0), (6) 

where i,j, k are the color indices, a and (3 are spinor indices, and P is the momentum 
of the diquark. Note that the flavors of the quark fields are not written out explicitly, 
so the flavors of the quarks could be the same or different. Since the diquark (like an 
anti-quark) must furnish the representation 3 of the color group SU{3)c, we can define 
a colorless BS wave function by Xp{xi^X2tip = j\£^^^Xp{xi-,X2)ap , 

xM,x,U = 6^^''=(0|TC(^i)^^(x2)|P,A;) = e-^^y"^X,(pWe-^% (7) 
X,(x2,xi)^, = e'■^'{P,k\^i;'^ix2r^Pl,ix^r\0). (8) 

4 



where X is the coordinate of the mass-center, p and x are the relative momentum and 
the relative coordinate, respectively. As usual, we start from a four point function, 

Six„x,-y„yX^^^l,^ = (0|TC,(^i)C(^2)^^:(2/2)*<(yi)lO) . (9) 

The Fourier transform of this four point function is defined by (with indices suppressed) 

S{xi,x2;y2,yi] 



r d Pd P'd pd p' ^-^iPx+iP'Y-ipx+ip'y ^ p p/^ 
J (27r)^^ 



^^^^^e-^^^'^-^y^^^^^^'y Sp{p,p') , (10) 

where, ignoring the effects of the isospin violation, p = {pi—p2)/2, p' = {p^—p^)!^, P = 

P1+P2, P' = p[+P2, and pi, P2, p'l, P2 are the momenta of -ipaii^i)^ Vaii^'^)^ V'^U^i)*' 
^'^2(^2)*, respectively. The relative coordinates and coordinates of the mass-center are 
defined by x = Xi — X2, y = yi — y2, X = {xi + X2)/2, and Y = {yi + y2)/2. Then the 
BS wave function of a bound state composed of two quarks satisfies the following BS 
equation: 

x,(p)2S = / "^sPip, kf:£^^,Kp(k, ytj^iy^kti^^ , (11) 

where Kp{k,k') is the Fourier transform of the two-particle-irreducible kernel of the 
four point function ([9]) . In terms of the colorless BS wave function Xp{^iy^2)ai3 , this 
equation can be written as 

1 r d^kd^k' 



■f ■/ ■/ ■/ 



xApu^^ = 6 y mXiyM^'^m ' (12) 

where ^J)^ = ^Jl^J? - . 

J1J2 Jl J2 J2 J\ 

To obtain the normalization conditions for the BS wave functions we also need an 
inhomogeneous equation for the four point function, 

sp{p. T^r^l^,. = sp {p, p'r^^i^u (13) 



+ 



d'^kd^k' r.(o)/ i.\hi2,i2i[ (1 ufY'AJU'i q (u „/\iii2J2ji 



Near the diquark pole, we can isolate the contribution of the "bound state" (diquark) 
to the above four-point function, 

AT, 

5^(p,p')^S., = E^ %o_^^_,, x.(^)"i"X(poai 

+ terms regular at P° = Ep , (14) 
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where Ep = yP'^ + "^diquark 'on-shell' energy of the diquark. The inverse of 

Sp\p^p'), Ip{p,p'), is defined to satisfy the following equation!* : 



Sf\p,k)Ip{k,p') 



aia2,/32/3i 



= - (27r)^ \5\p - p')5^,p,5^,p,5'^^H'-^- ± 5\p + p')5^,p,5o.,^,5'^^-5'-^^\ , (15) 

where '+' is for the scalar diquark and ' — 'is for the axial- vector diquark. We also need 
an auxiliary quantity, 



Qpip,p') 



1 



Kp(k,p')±Kp(-k,p') 



Operating this quantity upon Xp gives the normalization condition for the BS wave 
function near the diquark pole P° = Ep (for more details see, e.g. Ref. [11]), 
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^{2n)^ ^p(p)g^|-^^(P'P') - \ Kp{p,p')±Kp{-p,p') 



Xp(p') 

(16) 



where the spinor indices are suppressed. Since the kernel in our approximation is 
independent of (see the following subsections) the normalization condition is reduced 

toQ, 



(17) 



Now, let us turn to the discussion about the kernel for the interaction between 
quarks. The sucess of the potential model for mesons tells us that the strong interaction 
between a quark and an anti-quark can be modeled by two important interactions, a 
one-gluon-exchange part (with an effective strong coupling) and a linear confinement 



* In general, the inversion of Sp'^ipjp') will be more complicated. However, in our case (omitting 

the isospin violation), the following is enough. 

^ Eq. (fT7|) is different from the expression in Lurie's book [11]. The reason is that we have used a 

different convention for the normalization of one-particle state: (p|q) = {27t)^6^{p — q) • In deriving 

this normalization condition, we have used the fact that Xpi~P)0a = ^Xp{p)af3 for an iso-scalar 

diquark and Xpi~P)f3a — Xp{p)ap foi' an iso- vector diquark which will be shown explicitly in the 

following subsections. 



part. We assume that this is also true for a two-quark system. Furthermore, the 
parameter in the kernel of the diquark system can be related to that of the meson 
system by the so-called one-half rule [12]: the kernel for the diquark is just one-half of 
that for the meson. 

2.1 BS equation for the scalar diquark 

In the following we will discuss the BS equation for a Lorentz-scalar and isospin-scalar 
diquark composed of d and u quarks. The definition of the BS wave function of the 
scalar diquark can be written as 

From this definition we can see that % (— = —Xp {p)ai3- The free two-particle 
propagators reads 

+mji,p')^iP2,j2,-p'))}, (19) 

where S{pi) and S{p2) are quark propagators. The kernel arising from the one-gluon- 
exchange diagram is (we don't consider the effects of the isospin violation throughout 
this paper) 

^P^ yP^P )aia2,l32l3i — g 

X {{I'rTxlki^VTXii^^Ap' - p) + m,h,p') - (/92, j2, -p'))} , (20) 

where T"" are generators of the fundamental representation of the color group SU{3)c 
and A^jy is the propagator of the gluon field in Feynman gauge. The confinement part 
of the kernel is assumed to be 

^P^ \P^P )aia2Af^i 

= ^{(7°)LY^,(7°)3,i^(^'nP -pO + m,ji,p') - (/32, J2, -P'))} , (21) 

where (after imposing the covariant instantaneous approximation [9], for some expla- 
nation about this approximation, see the following text) 



where the second term is introduced to remove the infra-red singularity of the confining 
kernel near the points p'^ = and a small parameter /i is introduced to avoid the diver- 
gence in numerical calculations. To determine the constant c in the above confinement 
kernel, one can compare the non-relativistic approximation of the BS equation with 
the Schrodinger equation where the effective potential is used in the potential model 
for mesons. In the non-relativistic limit one can show that the effective potential from 
the BS equation is 

-ic 



V{r) 



2 as 



(23) 



Comparing Eq. ( !23l) with the effective potential in the meson case, V{r)^cson = + 
Kr, we can see that setting c = AiiTH is suitable from the viewpoint of the one-half 
rule [12]. 

With the interaction kernel in Eqs. (1201) and fl?Il) . the BS equation ffT^ becomes 
(using the relation Xp {-p)(3a = -Xp {p)ap) 



where /^(^s) 



XpSv) = [S{pi)^S{p,)]j ^2^^^ 

7^^ ® 7^ (j9 - p') + 1 ® 1 K^"'^ {P-P')]-Xp {p') , (24) 
^Qs (p-p'Y ■ "^"^^ convenience, we define a deformed BS wave function, 

Xp^ (P)aia2 = Xp^ {.P)ia2C~a^ = {CXp {.p))a^a2 , (25) 



2^2 



where C is the charge conjugation matrix. With this deformed BS wave function, the 
BS equation becomes 

d^p' 



Xpip) = [CS{p,)C-'^S{p2)] 



(2vr) 



C7^C-i ® 7^ K^'^\p -p') + l^l K^''\p - p')] ■ Xp {p') ■ (26) 



This equation can be written in a more usual matrix form 

d^ 



Xpip) = S{p2) 



{2n) 



YXp^ {p'fl, K^''^ + Xp^ iff K(^'^ S{-p^) , (27) 



where the superscript 'T' represents the transpose of the spinor index. In deriving the 
normalization condition of the BS wave function, we also need its conjugation defined 

by Xp ip) = C^^7°x (p)7°) which satisfies the following BS equation: 

</> <t> 

dy 



Xpip) = S{-p2) 



{2n) 



7" X„ ip')l, K^''^ + X. ip) K^"'^ Sip,) . (28) 
<j> <t> 
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From Eqs. (1271) and (!28l) . we can see that x_p {~P) and Xp (p)^ satisfy completely the 
same equation. Therefore, once we obtain the solution for x {p) we also obtain the 
solution for its conjugate by using x (p) = x_p {~p)'^ ■ 

Now, using the decomposition of a general matrix of Dirac fields 



1 



1 „.,.—,• 



+ 75(V^*75^^') - 757^(#^'757mV''')1 , (29) 

J f3a 

we can parametrize the BS wave functions as (note the fact that the intrinsic parities 
of the quarks are the same) 

Xp ipUa, = k/i + 757^(^;/2 + p;/3)/m^ + {-th,a,,P^p;^h/ml] , (30) 



where = p-Pcf, / is the longitudinal projection of p along the diquark momentum P^, 
p^M = _ p^p^jm^ is transverse to P^, (a = 1, . . . , 4) are Lorentz-scalar functions 
of J) 2 and p^, (T^i. = |[7m>7i^]- 

To simplify the BS equation ( !26l) . we impose the so-called covariant instantaneous 
approximation in the kernel [9]: = . In this approximation the projection of the 
momentum of each constituent in the diquark along the total momentum P^ is not 
changed. For a diquark at rest, this requires the 'energy' {p^ or pg) of constituent 
particle not to be changed. This approximation is appropriate since the energy ex- 
change between the constituents in the diquark is expected to be small when we use 
the constituent quark masses in the BS equation. Under this approximation, the ker- 
nel in the BS equation is reduced to Kp^ {p^ — p') which will be used in the following 
calculations. 

Furthermore, we define a three-momentum BS wave function as Xp {pj = / ^ Xp (p) 



- u - ^ 2-71 

the ^'-integration in Eq. (1261) can be carried out with the following redefinitions: 



j ^faip), a = l,...,4. (31) 

As usual, these four functions are not independent of each other. Their relations can 
be obtained by projecting the BS wave function with (see, e.g. Refs. [13,14]) 

At, = ^i^^, i/, = ^« +.".), = (32) 



where the energy is defined by £'1^2 = \J ~Pt^ + '^12- The square of the covariant 
'Hamiltonian' Hi 2 is the square of the energy, = -^12' the projection operators 
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satisfy A^A^ = , A^ AJ = , / = 1,2. Furthermore, we will take the quark 
propagators to have the form of the free one which can be written as 



—I 



+p^- Ei+ ie){m^/2 + p^ + Ei - ie) 
-m^/2 + p^ - H2 



■ {-m^/2 + p^+ E2-ie){-m^/2 + p^- E2 + ie)' 

Operating the projection operators on both sides of Eqs. (!33l) and (jMl) leads to 

1 



(33) 
(34) 



and 



m^/2 + Pi, - Ei + ie 
1 

m^/2 + p^ + El-it 
1 



Ar 



(35) 



^t^S{p,) 



-m^/2 +p^+ E2-ie 
1 



A+ 



A;; 



(36) 



-m^/2 +p^ - E2 + ie ^ 
Now, by multiplying A^^^ and A2^(f, on both sides of equation ( l26l) and integrating 
out the longitudinal momentum p^ along proper contour (s), we can obtain the following 
constraint equations: 

(A+f^)«/„^(A2f<^)a^a2Xp^(pj7a2C7ai =0, (37) 
(Arf0)«'i«i {^tfc^)a'^a2Xp^ {pJ^a^C^a^ = . (38) 

Substituting the parametrization Eq. ( !30l) into Eqs. ( !37j) and ( !38l) we obtain the 
following constraint relations 0: 



/3 = 



m 



f2 



(39) 



where we have defined m = mi = m2 and consequently uj = Ei = E2 . 

In addition to the above constraint relations, we can also obtain other two equations 
by operating Af /"^^ and Af /"^ upon both sides of equation 



Afl±S{pi)0Af^S{p2) 



i2n] 



The above two equations can be written in matrix form: A2 fd,Xp ip)'^C{-^i f4>)'^C ^ = and 
^t-fipXp (p)'^ C {-^1 f'p)'^ — . These two equations are in fact Unear combinations of independent 
matrices: 1, f^, |( , and ^0^^. • From the equations of the coefficients of these matrices we can obtain 
the following consistent solutions. 
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which, after taking the covariant instantaneous approximation and completing the 
integration J ^ on both sides, leads to (using x = C~^x ^iid written in terms of the 
matrix form) 



(41) 



(27r) 



Multiplying both sides of Eq. ( 1411) by C ^75 from the right and taking the trace over 
the spinor indices gives 



/i(pj 



cj(4cj2 - ml) J (27r)3 



m 



(m^+p ■ p'W^""'^ - 2mV(^s)l f2(p') \ , (42) 



d^p' 



^|-mm^(F(^f)+4\/(is))j^(p^') 
(m^ + p^ . p;)\/("^) - 2mV(^s)j 72(p;) 1^ . (43) 



where V^'^^'> = —iK^'^^^ and = _^f^(ig)_ At this point, let us define three-vectors 
and which satisfy p'^ = —p'^^ , p'^ = —p^ and p^' ■ p^ = —p^ ■ p^ §. Using these 
definitions, after completing the azimuthal integration, we can rewrite Eqs. fH21) and 
in the one-dimensional integration form, 

1 



/a(lpj) 



^ " b=l,2 

where a,b = 1,2, the 'matrices' are defined by 



d\p:\Yi ka^(ipj>ip:i)A(ip:i)-^.^(ipj>ip:i)A(ipj) ,(44) 



An = 2uj\Lo + AGo), 



A12 = -^(m^Lo + Li - 2m2Go) , 
m 



and the counter terms are 



Dn = 2uj^Lo 



D 



12 



2 



m 



m - IpJ )Lo 



21 



mm^Lo , L)22 = 2(m^ - |pJ^)Lo 



^ Since p^' and are four-vectors perpendicular to the total momentum , which is a time-like 
four- vector, we can always accomplish this. 
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where the definitions of Lq , Li, and Gq are given in Appendix El- The parameter /i 
will be removed in the end of the calculation by letting /i — > (in fact, taking fi to be 
sufficiently small is enough for pratical calculations). 



2.2 Normalization condition for the BS wave function of the 
scalar diquark 

Now, we will discuss the normalization condition for the BS wave function of the scalar 
diquark 0. From Eqs. ffT^ and ffTUl) we can define the inverse of S^^^ as 



IpiP, = 4^"^^^"^^ i^^mp - p') [S{p,)^X,, [^(P2)7Ti, • 

From Eq. (|T7I) we have the normahzation condition for the BS wave function, 

24 J (^^p/^)^2'3i^|[5(pi)7°];'/3j5(P2)7°];iJXp/p)aia, = 1. 

This equation can be recast in the matrix form, 
1 f d^p 



(45) 



(46) 



48 J (277)4 



-^Aci'xApfi'aipfsi-p,)-' 

I 4> <p 



-Ci'x,Xpfi's{p,)-'x,Xpf 



(47) 



where ^=(1,0). Furthermore, we can separate out the longitudinal momentum by 
using Eq. ([27]): Xp {pV = '^S{p2)Xp {pJS{-pi) and define 



n Pj^p" 

XpM) = 75 <! hipj + -r^M) - ^^/^^-32^^4(pJ \ , 



(4J 



where 



h2{pX 
hiip,] 



d^p^ 
(27r)3 
d^p'^ 
(2^ 
d^p'^ 
(2^ 



■ _ 2V'-'^\p^ - p[) + V'^^'Xp^ - p[)\ Up[) 

v^^'\p, - pf-^w:) ■ 



(49) 
(50) 
(51) 



On the other hand, from Eq. (|28l) and the discussion following that equation, we have 

Cl'Xp^ipfl' = C[x_p^i-pf]^C-' ^ ^5(-pOt^(pJ^^^5(p2) , (52) 
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where 

= l^'^P^^ + ^^^^^t) + ^^M^^^^(n)|75 • (53) 
With these definitions, we can write the normahzation condition as 
"S/ (^Tv{s(-p,)i/P.)'"''Sfe)?5fe)?,,(p.) 

-S(-Pi)^,,(p.)<"*S(p2)J,^(p,)S(-pi)i?| = 1. (64) 

After integrating out the longitudinal momentum and carrying out the trace calcu- 
lation, we have the following one-dimensional integration equation: 

2m^co'^(/ii)^ 4mm0|pJ^/i2/i4 - ?7im<^(m| 4a;^)/ii^2| = 1, (55) 
where E^j, — P^j, ■ ^ — \s the energy of the scalar diquark. 

2.3 BS equation for the gixial- vector diquark 

Now we will derive the BS equation for the axial-vector diquark. Since we will not 
take the isospin violation into account, we can take the diquark composed of uu with 
— as example in the following discussion. The BS wave function is defined by 

X^;^^{xi.x,U = 6^^^0|T<(xOK^^(x2)|P^,r,^) = e-^-^| ^X,^ (p)a/3 6"'^% (56) 
Xt\x2,x,)p^ = e'^'{P^,r,k\Tu^{x2rui,{x,r\0) , (57) 



where i, j, k are color indices, r is the index of the polarization vector of the axial- vector 
diquark. The free two-particle propagator reads 



„0li2j2 
J 02/92 



<Hp,p')L\ri^,, - (27r)^{5^(p-p')[^(Pi)/]Sj^(P2)7° 

-Wi^P2,ji^j2,p'^-p')}. (58) 

The kernel arising from the one-gluon exchange diagram is 

z?(is)('„ „/yi«2,i2ji _ (^5's)^ 

X {(7V)L\^(/7^)aA,.(p' -p) - (A ^ P,J, - j2,p' - -P')} . (59) 
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On the other hand, the confinement kernel is assumed to be 



From the above kernel and S^^^ we can obtain the BS equation for the axial-vector 
diquark, 



{2n) 



C-f^'C-' ® 7^ K^'^\p -p') + l®l K^''\p - p')] ■ xt}{p') > (61) 

which has the same form as that for the scalar diquark, Eq. fl26l) . Similar to the case 
of the scalar diquark, we can parametrize the BS wave function of the axial-vector 
diquark by several components (?j (/ = 1, . . . , 8), which are functions of p^ and P<^: 



xtip) 



pppf'^+pPpf'^ 

t |y^2 t t ^2 J 



pPp^'P"^ + gP^'P''^ + gPPp''-^ 



a 



(62) 



where ep"* is the r-th polarization state of the axial- vector diquark, which satisfies 

Sr) pp 



s^p^PP = . Performing the same procedure as that for the scalar diquark in the previ- 



ous subsection, we have the following constraint relations for the coefficient functions: 

~ ~ „ ~ P^93 + mlgr ^ m ^ 
92 = 9& = ^, 9i = im ' 96 = —95- (63) 



With these constraint relations, the BS equation for the axial-vector diquark is com- 
posed of the following four integral equations: 



1 



97 



2u;m'^p^^{m^ - Au^) J (2 



d^p' 



mm,, 



2/2 / l\2 

K p, - (k-pJ 



9i 



95 



2uj' 



'pX-{K-pX]iy^''^ + ^y^''^)93 



(64) 
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^3 



^7 



fl^3^ (65) 



94 



2m 



l\2 



2mujp^^{m?^ - 4cj2) J (27r)3 



y{cf) 



^4 



^5 



- m^p^p^'^ ^ 2py[^{p^ ■ p[) + 3m2(p^ ■ p^')^) ^^^^^ 
+ 2 - m - Apy[^{p^ ■ p[) + 3m2(^» ■ p'f^ ^^^^^ 5^3 

irt t 



(66) 



p>f -(p^-p;)^l(\/(^^) + 2i^(^s))^3 , 



^5 



(67) 



where uj^ = m? —p^^= m? + |pj^ • After carrying out the azimuthal integration of the 



three-momentum p^', we have 
-1 



9a 



2u;|p|4(m2 -4cu2) 



d\^[\ 5a.(ipj,ip:i)^,(ip:i)-a,(ipj,|p;i)^,(|pj) ,(68) 



where _Bafe(|p I, Ip'l) and Ca^dp I, Ip'l) (a, 6 = 3,4,5,7) are defined in Appendix Rl. 
Using the Gaussian quadrature method to discretize the integral equations, we obtain 
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the following linear equations: 

97 = U^g^ + U^g^ + U^g^ + Uygj, 
= R393 + ^4 94 + R5 95 + R797, 

= T-iQ^ + T^g^ + T^g^ + Tjg^ , (69) 

where f/j, -Rj, Si, Ti (z = 1, . . . , 4) can be read off from Eqs. fl6^ . fl65|) . fl66l) . and fl67j) . 
respectively. After some algebras we have the following eigenvalue equation for Qj-. 

97 = Ut g-j — Us [r45"'^T43 — S^^Sas] [T^^^^T^j — 5*45^ 5*47] g^ 

~ ^4 [^53^^54 - 'S'53^S'54] [T'53^r57 - 5*53^ 5*57] g^ 

~ U5 [T34^T35 — 5*34^ 535] [T34"'^T37 — S'34"'^S'37] g^ , (70) 



where, for convenience, we have defined 

Tij = Tr^Tj - Ri^Rj , 5*^^ = S^'^Sj - R^^Rj . (71) 
After gj is solved out from Eq. (!7(7|) . we can obtain ^34,5 by the following equations 

93 = — [^45^^43 — 5'45^ 5*43] [r45^T47 — S'45^S'47] , (72) 
9i = - [^53^^54 - 5*53^ 5*54] [T^^^T^j - S^^^S^j]]]^ , (73) 
g, = -[T-,'Ts,-S^,'Ss,]-\T,-,'nr-S^,'Ss7]gr. (74) 

2.4 Normalization condition for the BS wave function of the 
aixial-vector diquark 

The inversion of the 'free' four-point function (l58l) can be defined to be 

Ii^iP,P'r:^:^%. = ^^"^■^'^"^■^(2vr)^5^(j> -p')[5(j>i)/].-^,J5(j>2)7^i, , (75) 

From this equation and using the fact X^Hp)a0 = X^\^p)0a we arrive at the normal- 
ization condition for the BS wave function (see Eq. (ITTl) ). 



-Cj'xl:\pf^'S{p2)-'xl:\pf^ } = l. (76) 
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Similar to the case for the scalar diquark, we define 



where 



-^g^^p;ila,^-^pXP;^a,J., 



m. 



m 



1-1%' -p:a^ 



• ■I'tf •■•'ip 

and gi]^', . . . , 5*7' are defined by the following integrations: 



_ If dy 3{p^ -p'Y -pfp' 



p'^p^ — (p ■ p'Y 



(cf) 



2m^ 



p'2 2 _ / l\2 



Then the normalization condition can be written as 

^,Trl^S{-p,)^^;^{pj(^^S{p,)^S{p,)^ 



(r) 



ip) 



-S{-pi)xt!iPj^''^S{p,)xt!{pJS{-p,)^\ = 1, 



(77) 
(78) 



(79) 



(80) 



(81) 
(82) 
(83) 
(84) 
(85) 
(86) 



(87) 

which is similar to the case for the scalar diquark. After integrating out the longi- 
tudinal momentum p^ and carrying out the trace calculation, we have the following 
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one-dimensional integral equation: 

2E^ r \pfd\^ 



f iPt r%J [. 



+m^|pj {m^ + 4:u; )g/g4'-m\pj {m^ + Au; )g^'g4' - Amm^\pJ g^%' 
-4m^ml\pJ%%' + mml\pf{ml + Auj^)g^'g^ + 2mJ|pj2(mJ + Auj'^)g^g^' 
+4mm^^\pfg^%' - Zmm%{m% + 4u:'^)gig'l - m%\p;^{m% + 4u:'^)glg^ 

+mmj|pj2(mj + ^uj"^)])^])^ - Smml\pfg^gQ - 4mluj'^\pfgQ%' 
where — • ^ — is the energy of the axial- vector diquark cp. 



3 The effective interaction of diquarks and the pion 

As pointed out in the introduction, to study the decays of baryons in the diquark 
picture, wc should first calculate the matrix element (7r(g)|T0*(x)(^*(?/)|O) . This can 
be determined by the effective low-energy interaction of diquarks and the pion. We 
will calculate the effective coupling constant in this section by presenting the decay 
amplitude of a process, where the axial-vector diquark decays into the scalar diquark 
and a very soft pion, in terms of BS wave functions of diquarks and (p . 

Let us first introduce the effective interaction vertex of diquarks and the pseudo- 
Goldstone-boson pion: 

= G^^^ J2 <l>'d'^' + h-c. , (89) 

where i is the index of the (anti-)fundamental representation of the color group SU (3)c, 
'K^{x) is TT meson field with h {— ±, 0) being the isospin index. Then the decay amplitude 
can be calculated out to the lowest order, 

— 3G 

(0^(P,)7r«(g)|(^"'^(P^), r) = (27r)^5^(g + - P^)— ^-=^^ ^'^^^(P^) , (90) 

where the summation over the repeated color index i is assumed and no summation is 
assumed for the repeated isospin index a. On the other hand, we can present this decay 
amplitude in terms of the following transition amplitude with the aid of the partial 
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conservation of the axial currents (PCAC) [15], d^A"^{x) ~ -\/2m^/^7r"( 



X) 



'-{P^)\Al{q)^\^'^'\P,\r), (91) 



where /^r ~ 93.3 MeV is the pion decay constant, and the axial-vector current is 
^'iii^) — "4^(11^1 iJ.'^u{z) ■ Furthermore, momentum conservation and Lorentz invariance 
lead to 

q^{<P\P^)\A-{q)\^^\P^),r)^2E^^ = P^ ■ (2vr)V(g + P^ - P^) A{q^) , (92) 

where q = P^ — P^p and • e^^^ = have been used. The factors ■\/2E^ and ■\/2E^ 
are introduced in Eq. ( l92l) for latter convenience 0. This can been seen from the 
normalization equations ( l55l) and (1881) : with these energy factors absorbed into the 
BS wave functions of the diquarks, we shall only normalize -yj 2E^ Xp ^ind \/2E^p Xp\ 
which are Lorentz covariant quantities, instead of y and y , respectively. 
Combining Eqs. (EOj), (EI]), and ([921) gives 

This equation shows that the amplitude A must develop a pole at = . As usual, 
we can decompose this amplitude into two terms: one has the desired pole, the other 
is regular: 

A{q') = B,iq') + B^iq') , (94) 

where Bi and B2 are smooth functions without extra poles at q^ = , and their 
dependence on q^ are expected to be very weak and hence are nearly constants when 
q^ is small enough, e.g. q^ G (0, m^) . With these expectations in mind, we have the 
following Goldberger-Treiman-like relations: 

B^{q^) ^ 3V2UG^^^ , B2{q^) ^ -3V2f^G^^^ . (95) 



Then, the effective coupling constant can be calculated approximately by 

1 



^ In fact, by including these energy factors in the definition of the transition amphtude we make A 
be a Lorentz invariant quantity since we have used the normahzation convention for the single particle 
state, (p|p') = (2^)2<53(p_p'). 
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The transition amplitude A{q'^) at = is regular and it is this value that will be 
calculated in the following. 

To calculate A , one can choose a specific polarization vector. Choosing e^^^ = 
and using Eq. (lU^ . we have 

(27r)V(g + P^-P^)^(g^) = ^^2^^. (97) 

Let us now turn to the calculation of the matrix element {(f)''{P^)\A~{q)\ip'^'\P^),r) , 
which can be represented in terms of the BS wave functions of the diquarks and ip 
(using Eqs. (152|) and ( ITTj) ). 

{<P\P^)\A-{q)W^^P,).r) = -{2^Y5\q + P^ - P,) j -^{27^Y5\p, ~ p',) 

X Tr |5(-p'i)J,/p:)^'^5(p2)?^^\pJ5(-Pi)7;.75 

= ^{27^f5\q + P,-P,) j ^Tr|5(-p'i)t/p:)^^^5(p2)?^;;\pJ5(-p07M75}, 

(98) 

where pi = P^/2 + p, P2 = P^/2 - p, p'l = P(j,/2 + p' , p'2 = P(p/2 - p' , and the delta 
function 5^(p2~P2) imposes the following constraints upon the variables in the integral: 

, _ -IpJIP^I cosg + ml-P^- P^+p.Pj 

= -\vf-\p^\Pl\cosd-\Pl\^/A+{p+Pl/2-mj2f-p[\ (100) 

p[-Pl = -\p^\\Pl\cose- {1/2 -p'Jm^)\Pl\\ (101) 

P[-P, = -|pJ'-(l/2-p;/m^)|pJ|P*|cos^, (102) 

p[-P^ = p^m^- ml/2 + {1/2 -p'Jm^)P^-P^, (103) 

where Pt=p- P^/m^ , p'^ = p' ■ P^/m^, P, = P - ^P^p 5 P[ = p' - ^P<P^ and cos 9 is the 
azimuthal angle between p^ and P^ . By momentum conservation, jP^I and P^ ■ P^ can 
be written as: 

P,-P, = "'^f~'' ^Pim,, \n\' = -ml + {P^)\ (104) 
where we have defined P^ = P^ — :^P^ . 
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To integrate out the longitudinal momentum p^, we decompose the propagators 
S{-p[) , S{p2), and S{-pi) into 

(/i - m) 



S{-P[) 

S{P2) 

S{-Pi) 



—I 



2ijjq [p^ + Pf- - ujq + ie p^ + - + Ug - ie 



Pe - + ujp~ie 



P, +m^/2 -ujp + ie 



AtipJ + 



p^ - m^/2 -ujp + ie 



p^ + m^/2 + ujp-ie 



,(105) 
(106) 
(107) 



where p[ = p^ + ^P, + Pl+{^- \)P, , u:. 



+ |pJ2+|pt|2 + 2|pJ|Pt|cOS^^, 



+ I pj 2 , m is the mass of the constituent quark within the diquarks and . 

~ P ~ P 

The modified projection operators are defined by = and Af = C^^(A^)"'"C^ 

and can be written out explicitly as 



2m„ 



± 



2ljJn. 



A?(pJ 



(lOJ 



v^p ^^p 

Carrying out the integration over the longitudinal momentum p^ , the transition am- 
plitude (l98l) becomes 



\P^)\A-{q)\^^\P^),r) = -{2nr6\q + P^~P^ 



d^p^ 



Tr(++1 



X 



\p^=m^/2-ojp 



Tr(+- 



+ 



(27r)3 

Pl,=m^/2+u}q-P^ 



{2Up - m^) (P| -UJp- tOq) {Up + UJq- Pi) (m^ + Up + Uq- P^) 



Tr 



\p^ = -m^/2-LUp 



Tr(++) 



+ 



\p^=—m^/2+ujp 



{2ujp + m^) {m^ + ujp + ujq- P^) {2ujp - m^) (m^ -Up-Uq- P^) 



Tr I 



\p^=m.^/2- 



Tr 



"1-^ 4, 



\p^=mip /2+u}p 



{Up + ujq + P^) (m^ -Up-Uq- P^) {Up + ujq + P^) (m^ + 2ujp) 
where Tr (±±) are abbreviations for the following expressions: 



(109) 



Tr f±±) 



Tr 



2uJr, 



(110) 



The traces are too lengthy to be expressed here (and so are the transition amplitudes). 
Therefore, we will give only the numerical results in the following. 



Numerical results 



21 



Motivated by the studies of mesons in the potential model [16] and in the BS formalism 
[9], we take the parameter k = 0.2 GeV^ in the effective potential ( 123|) and vary the 
effective coupling constant to find the solutions of the BS equations. For each value 
of the diquark mass, takes a certain value. On the other hand, from the analysis of 
the spectrum of the heavy baryons, Ref. [17] shows that the mass of the scalar diquark 
is related to that of the axial-vector diquark by — ^ 0.210 GeV. We find that 
for the scalar diquark, (m<^(GeV), a,) = (0.70,0.590), (0.75,0.570), (0.80,0.555), while 
for the axial-vector diquark, (m^(GeV),a,) = (0.91,0.490), (0.96,0.307), (1.01,0.093) . 
With the solutions of the BS equations as the input, the transition amplitude in Eq. 
(11091) can be calculated out. Then, from Eqs. ( l96l) and ( 1971) we obtain the values of 
^(0) and Gt,^^p . The results are listed in Table [H It should be noted that the sign 

Table 1: The transition amplitude A{q^) at = and the effective coupling constant 
Gtt^ip corresponding to various values of the mass of the scalar diquark (the mass of 
the axial- vector diquark is related to that of the scalar diquark by m^p—m^ ~ 0.21 GeV). 



(GeV) 


0.70 


0.75 0.80 


^(0) (GeV) 


1.32 


1.43 -1.53 




3.35 


3.62 -3.88 



of the coupling constant Cj^^^p is an artifact in numerical calculation (which is related 
to the arbitrariness of the sign for the BS wave functions). A physically observable 
quantity, e.g. the decay width discussed in the following, depends only on the norm of 
the effective coupling constant IGtt,^;^!. 

At this point, we want to point out that the numerical results show that the "++" 
component (the sum of the first and fourth terms in Eq. (11091) ) gives the most important 
contribution and the "H — " and " — h" components give less important contribution 

(< 15%) to the total transition amplitude. The " " component (the sum of the 

third and sixth terms in Eq. (11091) ) gives very small contribution to the total diquark- 
transition amplitude. 
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4 Strong decays of heavy baryons 



In this section, we will turn to the calculation of the decay widths of the processes 
TiQ^ Aq + tt . As pointed out in the introduction, the heavy baryons Ag and Sg^ are 
regarded as bound states of the heavy quark Q and the diquarks, which are composed 
of two light quarks. In this picture, we can express the decay amplitudes in terms of 
the BS wave functions of Aq and Sg^ (taking Eg^ as an example), 

<A$(i^A)vr+(?)|4^(i^E)> 

d'^{xiX2yiy2Uv)Xp {x2,Xi)Sq{Xi - yi^^Xp ,a(i/i,Z/2) 

X A-/ix2 - u)A^''-^\v - y^) 5^(vr+(g)|T0X«)^:M|O) , (HI) 

i 

where the superscript i is the color index, 0* is the field of the scalar diquark and is 
the field of the axial-vector diquark. A_^^ and A"^''^'^ are the inverse of the propagators 
of the scalar diquark and the axial-vector diquark, respectively. They are defined by 
i(x, y) A^(y, z) = 5\x - z) and A-^'^^x, y) A^,A^(y, z) = 5^^5\x -z). Sq{Q = c, h) 
is the quark propagator. The BS wave function of the heavy baryons are defined by 

X, (X1,X2) = (O|TV^^(xi)0^(x2)|Pa), (112) 

X,^ ..(2/1, 2/2) = (0|T^^(yi)<(y2)|PE). (113) 

If not explicitly pointed out, the summation is understood for repeated color index i . 
The complex conjugate of the BS wave function is defined by Xi) = x (a^i, 3^2)*7° • 
To simplify the analysis, we take, as usual, the propagators of the diquarks and the 
heavy quark to have the forms of the free ones. Then we have 

A^\x,y) = iiD, + ml)5\x-y), (114) 
A~''''\x,y) = t[{n, + ml)g''' + d:d^]s\x-y), (115) 

SQ(x,y) = [ 7^- e-'f", (116) 

'^^ J (27r)4 ^-mQ + te ' ^ ' 

where mQ is the constituent mass of the heavy quark which can be determined by, e.g. 

fitting experimental data to the results of the potential model for mesons, and m,:^ 

are the masses of diquarks and (p , respectively. 

The matrix element (7r+(g) |T (j)^{u)Tp^{v)\0) can be calculated out to the lowest order 

with the effective interaction vertex defined in Eq. 0891) . 

{7r^{q)\T<P\umvm = -3G^^^-^ J d'z e^"^ A^{u - z)A^{z - v),, , (117) 
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where the factor 3 comes from the summation over the color index. Then Eq. fillip 
becomes 

—3G 

(A$(PA)7r+(g)|S++(Ps)) = -^(27r)V(g + Pa - P^) 

(2vr)V(p; -pi) Sq{p,)-%^Ap) , (118) 

where pi is the momentum of the heavy quark Q within the baryon Sq and Pi is the 
momentum of the heavy quark Q within the baryon Ag, p2 and are the momenta of 
the diquarks (f and , respectively. These momenta are related to the total momenta 
of the bound states by the following equations: 

p[ = X[P^+p', p'^ = X'^P^-p', Pi = AiPe+P, P2 = X2P^-P, (119) 

where the parameters Ai^2 and A'^ 2 are defined by 

A, = ^^. A. = ^^. (120) 

Following Refs. [6] and [7], we can present the BS wave functions of the baryons as 

X, (P) = iu^{v')N{p[,p'^)A^{-p',)Sq{p\) , (121) 

A t t 

where ma is the Dirac spinor of Ag . Since there is a delta-function 5{p'^ — pi) in Eq. 
(11181) . to the leading order in l/mg expansion, we have 

X, (pO Sq{p,)-' = -— ' u^{v')N{p[,p[) . (122) 

A [p^ + ms - /a) - + * 

Similarly, to the leading order in l/mq expansion, we can present the BS wave function 
in the baryon Sg^ as 

<W ^ (123) 

where m = 1 is for Eg , m = 2 for Eg , Eq is the binding energy of the baryon Eg-* 
which will be determined in Sect. El For more details on discussions about the BS wave 
functions of heavy baryons, the readers are referred to the original references [6,7]. For 
convenience, we also give some results in Appendix [Bl In the above equations, we 



have defined the following quantities: Wg = y ~(Pt ~ Pa)"^ + > = y~P^ + > 
P'l = p' ■ v' — X^rriA , Pi = p ■ V — \2'm-^ . p'^ = p' — {p' ■ v')v' and p^ = P — {p ■ v)v 
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are perpendicular to the "velocities" of the baryons, v' = P/^/m^ and v = Pj^/nij^, 
respectively. = {m\ + m| — q^)/ (2mE) = m^t> ■ t>' is a constant, = — is 
perpendicular to Ps • The polarization vectors of the baryons Eg and Eg are given by 



B^\v) = ^{^, + v,)^,u^{v), Bf\v)=u,{v), (124) 

where Uy^{v) is the Dirac spinor and 'U^(f) is the Rarita-Schwinger vector spinor. 
B^^\v) satisfies the following conditions: 

fB^^\v)=B^^\v), v^Blr\v)=0, rBl^\v)=0. (125) 

The above constraints for m = 1 can be seen from ^Ms(f ) = uj:{v) while for m = 2, they 
are the properties of the Rarita-Schwinger vector spinor with spin |. The expressions 
of N{p',p') and M^''{p ,pj are given in Appendix [B] (these expressions are extracted 
from Eq. (12) in Ref. [6] and Eq. (28) in Ref. [7], some notations in this paper are 
different from those in the original references). On the other hand, the delta- function 
6{pi — p[) leads to the following relations: 



p^.pl = -m^lpjVn^ -IcosO, (126) 



p'^ = VLp^ - - icos^ + firrtE -mA, (127) 

= - IpJ 2 sin^ e - \vt\p^ I cos e - (ms + V^^^ - 1 1 \ (128) 

where Vt = v' ■ v . 

The poles of p^ and are all from the prefactors in Eqs. (11221) and (11231) since 
N{p'^,p'J depends linearly on and M^''{p^,pJ contains up to second order [6,7]. 
After carrying out the integration over p^ in Eq. (I118p . we have 



where 



rP(n^)= /^/^f^ q,M-^{p}^{p'}\,-w^ 

^ ^ y (27r)3 \ 2W,{W, + + m^) [{W, + - P^f - W^] 

QxM''{p)N{p%^^_^^^^^^P^ 

^ 2Wg{-mj: + P/ + Wg + E^ + m^) [{-m^ + P/ + W^Y - W:fi 

[(-m^ - E^ + ms - Pi? - W^] [{m^ + E^f - %V ' 

25 



The tensor function M^^ is defined by 



M^p{p) = g^PM.ip) + -^M,{p) - -^Msip) , (131) 



where g'^'^ = diag{l, — 1, — 1, — 1} is the Lorentz metric tensor, M„ {n = 1,2,3) are 
defined in Eqs. (IB6P - flB8P in Appendix [Bl When calculating the contraction, we will 
use the orthogonal condition v'^Blr\v) = 0. Furthermore, the integrations of terms 
containing p^B^fT^ {v) have the following form: 

j d\ p^B^;^\v)hip^,Pl) = P'/B^;-\v)gi\Pl\) , (132) 
where h is arbitrary smooth function and 

9i\n\) = I d\^p^h{p^,Pl). (133) 

Then only terms containing the factor P^^'Bp"^^ (v) can appear in the final expression. 
The calculation is straightforward and we have 



(ms - mA^)m^\p^\ 


COs6M2a — 


Pl\{mlM,a + \p, 


pCOS^ OM^a) 


2Wp{Wp + 1 
(ms - mA^)m^\p^ 


^^ + m^)[iWp + m^-Piy- 
\coseM2b-\Pl\imlM,b+\p,\ 


2 cos2 OMsb) 


2Wg{-m^ + Pi + Wg + + m^) [(m^ -W^- 
{rriY. - mA^)m^\pJcoseM2c - |PaI("^^^1c + \pj 


Ply - wi\ 

2 COS^ ^Mse) 



where A^a = f^{p')\p^=W^ , A";, = A^(p')lR,=-ms+W,+P* ' = A^(p')lp,=-m^-BS , ^na = 



M„(p)|p^=vyp, = M,,{p%__^^^^^^p^, Mnc = M„{p)\p^^_^^_E^ , {n = 1,2,3). 
Furthermore, \Pl\ = m^V^^ — 1 and Pi = uia^ are constants. Note that the energy 
factors \plEl and \/2Ey, will be absorbed into the corresponding BS wave functions. 
This is convenient since they appear in the normalization conditions ( lC4p and ( ]C9I) 
given in Appendix O and make the normalization conditions be Lorentz invariant. 
In the rest frame of t}q , the differential decay width of tIq' — > Aq + vr is 
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where dfl is the sohd angle of the particle in the final state, and P^y is the three- 
momentum of the baryon Aq . The Lorentz-invariant amplitude in our case can be 
written as 

M = -3G,^^V{q') u^{v')B^;^\v)P^r (136) 

We will calculate the unpolarizcd decay width averaging over the spins of the initial 
states and summing over the spins of the final states. For Sg , we have 



\ \Mv\ s)Bf\v, 5)P;'f = ^ (fi - 1)(1] + If , (137) 

s',s 

while for Eg , we have 

\ J2 \Mv', s)Bf\v, .)Pi'f = ^ (0 - 1)(0 + If . (138) 

While deriving the above two equations we have used the following formula for the spin 
sum of the Dirac spinor: 

UB{kB, s)uB{kB, s) = ^'^^ ^ , B = AqOtT.q, (139) 

s=l,2 

with Vjy^^ = v', v^^ = V , and the formula for the spin sum of the Rarita-Schwinger 
spinor [11, 18], 

J2 s)u^{K s) = ^ [-9''' + + Iv'V' + ^(7^^^^^ - Yv^) \ , (140) 

s=l ^ ^ 



where the fact (^+1)^ = 'ji +1 has been used in the final step. The total unpolarized 
decay width in the final form then reads 



r(E« ^ + tt) = -|- m')? - l)(n + If . (141) 

5 Numerical analysis 

We will first discuss the parameters appearing in this paper. These parameters include 
m<^, rriip, Kb, mg, and Eq . The study of mesons in the BS equation approach [9] shows 
that the values of the masses of the heavy quarks rric — 1.58 GeV and = 5.02 GeV 
lead to predictions in good agreement with experiment. 
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Now, we discuss the ranges of the masses of diquarks and (p. In contrast with 
the colorless hadronic states, diquarks (and other color states) are not free. Therefore, 
their (effective) masses can not be measured in experiments. In our discussion, we will 
treat these masses as parameters varying in reasonable ranges which satisfy various 
constraints from physical considerations. From the analysis of the spectrum of the 
heavy baryons, Ref. [17] shows that the mass of the scalar diquark is related to that of 
the axial-vector diquark: — ^ 0.210 GeV 0. In our calculation, the mass of 
is taken to be in the range G (0.70, 0.80) GeV and the corresponding range of the 
mass of ip is G (0.91, 1.01) GeV. 

Now consider k,^ . It is argued in Refs. [6, 7] that the parameter k,^ in the effective 
potential (See Eq. flB3p ) can be ranged approximately from 0.02 GeV^ to 0.1 GeV^ . 
Furthermore, by studying the average momentum of 6-quark in A;, and comparing 
with the value of this quantity derived from the experimental value of the average 
momentum of the 6-quark in the B meson with the aid of HQET, the authors in 
Ref. [10] show that can be constrained to a narrower range: "When are 0.7 GeV 
and 0.8 GeV, are roughly in the ranges (0.02-0.06) GeV^ and (0.02-0.04) GeV^ 
respectively." Following this, in this paper, we will calculate the decay widths in the 
range G (0.02-0.06) GeV^ 

Now we will determine + Eq. In the heavy quark limit, the baryons Sg and Sg 
should be degenerate and the dynamics inside them are the same. In the heavy quark 
limit, we can write out the masses of the baryons as 

--niQ + m^ + Eo + O (^^^ , (142) 

where m^ + Eo is the value to the leading order in 1/mg expansion and then is universal 
for all heavy baryons (with one heavy quark). Since mb 3> m-c, the l/mg corrections 
for 6-baryons are much smaller than those for c-baryons, hence we shall take the data 
of 6-baryons as the input to calculate this quantity. From the recent results of CDF 
Collaboration [4] , we have lI 

m^ + Eo^ 0.81 GeV (143) 



II — 7710 = 0.211 GcV for the c-baryons and to^ — — 0.208 GeV for the 5-baryons (the masses 
of 6-baryons are taken from Ref. [4]). 

** We calculate the following quantity by using the spin-averaged mass of Ej, , m-^^ = {2m^^ + 
4my. )/6, where TTij. = {m + + m )/2 and niy, = (m -t- m , )/2. Notice that all effects of the 
isospin symmetry violation are omitted. 
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for Sg when the I/trq corrections are omitted. For Ag, we have 

m^ + Eo^ 0.60 GeV (144) 

when the l/mq corrections are omitted. 

With the parameters determined above, the decay widths of the processes ^ — > 
Ac,6 + vr can be obtained. The results are shown in Table [21 From this table, one can 

Table 2: The decay widths r(S^*^^ Ac,fc + 7r) to the leading order in l/mq expansion. 
The violation of SU{2) symmetry is not taken into account. The unit of is GeV, 
the unit of F is MeV, the unit of is GeV^ . The mass of (f is related to that of by 
— m^^ 0.21 GeV . 



m0 


0.70 


0.75 


0.80 




0.02 0.04 0.06 


0.02 0.04 0.06 


0.02 0.04 0.06 


r(Sc) 
F(s:) 


6.61 4.83 3.75 
18.88 14.83 12.11 


4.91 3.85 3.15 
17.54 14.26 12.01 


3.95 3.26 2.77 
15.99 13.49 11.69 


r(s,) 
r(E*) 


13.45 10.20 8.10 
17.74 13.76 11.09 


11.16 8.88 7.34 
15.76 12.68 10.57 


9.47 7.88 6.73 
13.92 11.65 10.00 



see that the theoretical result for F(S*) is consistent with the experimental data [2], 
F'^^P(S*) ^ (15- 16) MeV. However, the theoretical value for F(Sc) is bigger than the 
experimental data [2], F™p(Sc) ~ 2.2 MeV. We attribute this to the l/mq corrections 
which are not taken into account in this paper. To look at this more transparently, 
one can estimate roughly the l/mq corrections as follows. If the corrections to the 
magnitudes of the BS wave functions are Aqcd/'^c = 0.16 (for Aqcd = 0.26 GeV), the 
corrections to the final results of the decay widths can be very large, (F + 5F) /F = 
1.84 to 0.49 for c-baryons. For 6-baryons, Aqco/'^b = 0.05, then we have (F + (5F)/F = 
1.22 to 0.81. Since the corrections for 6-baryons are much smaller than those for c- 
baryons we expect that the predictions for the decay widths of are far more precise 
than those for T,i*^ . 
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6 Discussions and conclusions 



In this paper, we have first studied the properties of two kinds of diquarks, the scalar 
diquark and the axial-vector diquark ip, in the BS formahsm. We have derived the 
BS equations for these two kinds of diquarks and studied all the BS equations under 
the covariant instantaneous approximation, which allows one to obtain the BS wave 
functions in a general coordinate system directly. With these BS wave functions of the 
diquarks, we have calculated the effective coupling constant among the diquarks and 
the pion, G-,,^^. We find that this effective coupling constant is IGtt,/,,^! G (3.35,3.88). 
With this effective coupling constant, we have calculated the decay widths of the heavy 
baryons t!q {Q = c, b) in the BS formalism in the heavy quark limit mg oo. 

There are two parameters m<^ (and m^) and in our model which can not 
be determined in this paper. Following the arguments in Ref. [10], we take k,^ G 
(0.02- 0.06) GeV^ For the diquark -masses we take G [0.70, 0.80] GeV and G 
[0.91, 1.01] GeV. With these ranges of parameters, we give the predictions for the decay 
widths of Sq'' Aq + tt (the effects of the isospin violation are not taken into account): 

r(Sj ^ (2.77-6.61) MeV, r(S*) ^ (11.69- 18.88) MeV , (145) 
r(S J ^ (6.73 - 13.45) MeV , r(S*) ^ (10.00 - 17.74) MeV . (146) 

From Eq. 01451) . we can see that the calculated value for r(Zl*) is consistent with 
the experimental results [2], r'^''P(S*) ^ (15- 16) MeV, but the calculated value for 
r(Sc) deviates from the experimental results [2], r'^^P(Sc) ~ 2.2 MeV. We attribute 
the deviation to the l/rric corrections which are not taken into account in this paper. 
However, we expect that the predictions for 6-baryons are much more precise than 
those for c-baryons since rrif, ^ m^. Furthermore, the above results show that the 
decay widths of Sf, — > Af, -|- vr will be much bigger than those of Sc — ^ Ac + vr. 

For comparison, let us quote the results in Ref. [19], where the decay widths of 
were calculated in the bag model (in units of MeV): 

r(S+)=4.35, r(S0) = 5.65, r(S,-) = 5.77, (147) 

r(s*+) = 8.50 , r(sf ) = 10.20 , r(s*") = 10.44 , (i48) 

which are comparable with our results in Eq. (11461) . 

In this paper, we omit all the l/mq corrections in the calculations. If the l/rriQ 
corrections are taken into account, we can take the value of r(Sc), which has been 
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measured precisely in experiments, as the input to constrain the parameters k^, m^^ 
and m<^ and make more precise predictions for r(S*) and r(s[*''). The full l/mq 
corrections include the 1 / itlq corrections to the kernel, to the heavy quark propagators, 
and to the BS wave functions. The l/mq corrections for Ag have been discussed in 
Ref. [20]. Since the study of these l/mq corrections is very complicated, this is beyond 
the scope of the present paper and will be discussed elsewhere. 

Finally, the effect of the isospin violation is not taken into account to make the 
presentation of the calculation more transparent. After taking into account the effect 
of the isospin symmetry violation one can obtain more information on the properties 
of heavy baryons (if including the strange quark s, one can also give predictions for Qq 
and Sq). This will be the subjects in the future work. 
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A Some definitions used in the previous sections 

In this appendix, we give definitions of functions used in Subsections 12.11 and 12. 3[ The 
sixteen 'matrices' Bij in Subsection 12.31 can be written as [2 

= -2u;'\p\'\pfLo + 4\p\'\pfL, + 2{2m' + u;')L, 

- 4cu2|pnpfGo - 16\p\^\p'\^G^ + A{2m'' + co'')G2 , (Al) 

= mm^{\p\^\p'\^Lo-3L2), (A2) 

S35 = -2ml\p\\L^-2G^), (A3) 

537 = -4ml\p\\-\p\'Lo + L,-2\p\'Go-AGi), (A4) 



In (and only in) this appendix, p and p' are in fact p^ and p^, respectively, which appear in Sect. 
[2j The change of the notation is to make the expressions more transparent. 



31 



B43 — 

B44 — 
B45 — 

B47 — 

B53 

B54 

B55 
B57 

B73 = 

B74 — 



rriu 



m 



m^\p\^\p'\^Lo + 2\p\^\p'\^Li + 3mU2 



-2m'\p\'\pfGo - S\p\'\pfG, + dm'G^] , 



-2 



Amm^lpl'^ 



|p|% + Li-2G'i], 



2m^^|pP 

'^"^'-(Li-4Gi), 



m 



= \p\' [\pW\'Lo -L, + 2\p\'\pfGo - 2G2] , 



2m\p\ 



2 r 



rriu 



-\p\^\pfLo + L2\, 
4|p|^pLo + Li-2G'i] , 



2ml\p\\Lo + 2Go), 



B75 
B77 



- \\pWLo -L, + 2\pWGo - 2gJ , 

L J 
-2|p|^[m'Lo + Li-2G'i] , 
4u;^\p\\Lo + 2Go), 



and the counter terms are 



C33 = 4|p|^(m^ - IpHLo , C34 = -2mm^\p\^Lo , 



C35 = 2m^\p\ Lo , G37 = 8m^|p| Lq , 

C43 = -^^^(m^ - IpHLo , C44 = 4a;^|p|% , 
m 

C45 = -8mmJp\^Lo , C47 = ^^Lq , 

m 



(A5) 
(A6) 
(A7) 

(A8) 

(A9) 
(AlO) 

(All) 

(A12) 

(A13) 

(AM) 

(A15) 
(A16) 



(A17) 



(A18) 



C53 = C54 = , C55 = 4|p|^(m2 - |p|2)Lo , C57 = -2m^|p|% , (A19) 

C73 = C74 = 0, Cv5 = -2b|^K-|pnLo, Cr7 = 4b|VLo, (A20) 
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where the functions Lj and Gi {i — 0, 1, 2) are defined by 

Attk 



J (27r)3 



)3 [(p_p/)2+^2] 

/" d^p' AttkIp ■ p') 
/■ (iV 4:7rK{p ■ p'Y 
d-V 2(?2 1 



(27r)3 3 {p-p'f + iJ? 
d'^p' 2gl p ■ p' 



(27r)3 3 (p-pO^ + A*^ 



with 



_ 2« 



(27r)3 3 ip-p')^ + ii'^ 

\pf 



2fi\P'\) ^ 


1 d\p' 




j d\p' 




JdW 


-j(\p'\) ^ 


jdW 


-ji\p'\) ^ 


j d\p' 


-j{\p\) = 


j d\p' 



IMIpUp'I)/(Ip'I), 

\LmAp'\)f{\p'\), 

|i^2(|p|,b1)/(b1), 

ICodpUp'D/dp'l), 

|Gi(|p|,|p'|)/(|p'|), 

\G,{\p\M)f{\p'\), 



TT (|p|2+ |p'|2 + ^2)2_4|p|2|p/|2 ' 

K \p'\ j m\p'\ i\p? + \pf + f^') 



{\p\ - \p'\r + 



47r bl 1 + + ;,2)2 _ 4|p|2|p,|2 + + |y|)2 + ^ 



(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 

(A27) 
, (A28) 



^2 = ^^MpIIp'I 



(|p|' + b' 


2 + /i2)2_2|p|2|p'|2 


(b 


2 + b' 


2 + //2)2_4|p|2|p/|2 



+(br+ b'r+Ai')iog 



Go 
Go 



as b'l 



log 



b'l)' + /^' 



-\p'\? + ^' 

(|p| + |y|)2 + ^2p 



svr bl "(bl + b1)' + /^'' 



(A29) 

(A30) 

(A31) 

\p'\f + l^' 



(1^1 + 1^1)2 + ^2 

(A32) 



where as — gU (47r) , p • p' = — bl b'l ^ ■ 
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B Some results in Refs. [6,7] 



In the previous sections, we have used the BS wave functions of the heavy baryons Aq 
and TiQ^ in Refs. [6,7] as the input to calculate the decay amplitudes. In this appendix, 
for readers' convenience, we write out some expressions (taking from Refs. [6,7], but 
in a different notation) used in this paper explicitly. For more details, we refer the 
readers to the original references. 

Following Ref. [6], one can write the BS wave function of Aq to the leading order 
in l/mg expansion in terms of the product of a scalar function N{p') and the spinor 
of Aq, 

X, ip') = iA4-p',)S{p[)N{p')uA{v') , (Bl) 

A 

where the scalar function N{p') is given by the following integration: 

/d^k ~ ~ ~ 

^(V^i + 2p;^)0pJA;O. (B2) 

0/3^ depends only on the norm of and has been solved numerically in Ref. [6]. The 
kernel in the above equation is defined by 



[(p; - A;t)2 + /x2]2 V I y^, J (27r)3 (/2 + ^2)2 ' 

3 [{p:-hy + ^l'][{p:-h)' + Ql]' ^ ' 

where Ql = 3.2 GeV^ , is a parameter which is taken to be small enough so that the 
numerical result is insensitive to this parameter. 

Following Ref. [7], one can write the BS wave functions of Tjq^ to the leading order 
in l/mq expansion in terms of the product of a tensor function M^'^{p) and the spinor 
of Eq\ -BjL™^ (m = 1,2, corresponding to Eq and Eq, respectively). 

The tensor function is given by 

M^'^(p) = g^^M^{p^, \p\) + ^M^fe, \p\) - ^M3(p„ \p\) , (B5) 



where 
Ml = 



(2^ 



2p2 



{V, + 2p,V,) - C - ' V, ) , (B6) 
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Mo 



(2^ 



A + D 



p2 



C 



/ 2 2 \Pt ' k^ 



(p ■ hf ~ ■ 

p2 



(B7) 



d^kt 



(27r)3 
+ C 



■ P,-K~ ml{3{p^.h)^-pfk^)+2p^{p^.hy 
R^^Vi + ^1 V2 



+ D 



t t 
2pf 



iV^ + 2p^V2)+p, 



p2 



(B8) 



The three wave functions, A, C, and Z) in Eqs. (IB6p - flB8p are scalar functions which 
depend only on the norm of kt and have been calculated numerically in Ref. [7] . 

It is worth pointing out that the following conventions have been used here, p^ = 
\pj'^ , p^ ■ kt = \p^\\kt \ , etc., which are different from our conventions used to discuss the 
BS equations of diquarks. 



C Normalization of the BS wave functions of the 
heavy baryons 

In this appendix, we will give the normalization conditions for the BS wave functions 
of Aq and Sg^ . The normalization conditions of the BS wave functions of Aq and Sg'* 
have been discussed in Refs. [6,7]. In order to calculate the decay amplitudes (other 
than the weak transition amplitude) of baryons in our case, it is necessary to obtain the 
BS wave function for each baryon separately. We start from a normalization equation 
which is similar to that of the diquark in Eq. f|T71) . Furthermore, in the heavy quark 
limit rriQ 00, the normalization conditions in this section should take the same form 
as those obtained through the normalization of the Isgur-Wise functions at the zero 
recoil point, ^{Q = 1) = 1 [6,7]. This can be checked with explicit derivation. 



Normalization ofXr, 
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The normahzation equation for the BS wave function of the baryon Aq is given by 

where s and s' are indices of the spin of the baryon and 

IpMp) = l{27rmp-p')SQ\p,)A-\p2). (C2) 

Then we have 

i f d^p _ d 



3 



After carrying out the integration over p^ , the normahzation equation becomes (mul- 
tiplying 5s,s' on both sides and summing over the spin indices) 

it. I St^ wMi^ {^^' - -;)-^ - -^<-« - - • 

which can be reduced further to a one-dimensional integral equation, 

Ea fd\pj \pfQ>i-2Wp4>2) (~ r .ji 

a / TTT u/ ( J7 ^ ^TTT^ i 01 + 2 1 - 2Ai - 2A2 + m,^ 02 [ = 1, 

6mA i 27r^ PVp(Eo + - IVp)^ I ^ ^ J 

(C4) 



where Wp = W + , (pi,2 are defined by 



~ r d^k ~ ~ ~ r d^k ~ ~ 

MpJ = y (^ - K) K^k,) , 02(n) = y - K) ^pSk.) . (C5) 

As pointed out in the beginning of this section, in the heavy quark limit, the normal- 
ization condition above should reduce to that obtained by the normalization of the 
Isgur-Wise function at the zero recoil point. This can be checked easily. In the heavy 
quark limit, itlq 00, we have A'^ = 1 and Ag = 0, then Eq. ( 104^ becomes 

fd\pj \pfi^,-2Wp^2r 
3mA J 27r2 2Wp{Eo + m^ - Wp)^ ' ^ ^ 

One can see that this is the same equation as that given in Ref. [6] at the zero recoil 
point = 1 (see Eq. (26) in Ref. [6], notice that Q is written as uj there) ul. 



W In fact there is an extra factor E\/ (3mA) in this equation when compared with Eq. (26) of 
Ref. [6]. The energy factor EA/m\ is required to make the BS wave function be a Lorentz scalar 
in our convention of one-particle states, (p|p') = {2tt)^S^{p — p') , which is different from that used 
in Ref. [6], (p|p') = {27r)^5^{p — ■p')Ep/m. Furthermore, the factor 1/3, which comes from the 
summation of the color indices, does not appear in Eq. (26) of Ref. [6]. However, the ignorance of 
this color factor does not affect the results when one calculate quantities like those in Ref. [6] . 
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Normalization of xt 



For the baryon Tiq^ , the normahzation equation of the BS wave function is 



(27r; 



8 ^P. 



d 



where 



Ip^{p,p'),. = TT {27:Y6\p~p')Sq\p{)^:^\p,),^. 



(C7) 



Multiplying dggi on both sides of Eq. flC7p and summing over s and s' , the normalization 
equation becomes (using the relations given in Appendix [Bl) 

d^P, \2{v ■ v) 

(27r)3 18m^^W^{Wp - Eq - m^f 



X 



3mjMiM{ + ml\pf{MiM'^ + MaM^ - M^M[) + bJ^MsM^ 
+ SmJWf (Mio - PFpMn) [Mio + {W^ - 2^b)M^^] 
+ bj' [w;'(2M3i + m^M22) [W^p(2Ab - W,,)M^^ - A^Mio] 

+ W^p'(Mio - ABMn)(2M3o + m^M^i) + m^Mao [(A^ - 2Wp)Mio + W^M^ 
+ \pf [w^iMsi + m^M22) [Wp{2/\B - W^p)M3i - A^Mso] 

+ WliM^o - AbMsi){M3o + m^M2i) + m^Mao [(A^ - 2Wp)M^o + W^M^i] 

(C9) 

where A^ = Eq + , Wp = ^JJpJ^~-\-m^ , v ■ rj = E-^/raY., M„ (n=l,2,3) have been 
defined in Eqs. (IB6])-(IB8]), and M„^(nm = 10,11,20,21,22,30,31) are defined by the 
expansion of M„'s , 

Mi = Mio+p,Mn, M2 = M2Q+p,M2i+pfM22, M3 = M^o + P^M^i . (CIO) 

M^{n = 1,2,3) in Eq. ^ are given by M[ = M^, = M20 and = -(M30 + 
m^M2i) — PiiM^i + m^pM22)- In the heavy quark limit, Eq. (1C9P reduces to 

\pf{v-ri) 



<m 

27r2 18m4 W^p(H^p - - m^)2 

- 3mjMiM( + ml\pf{M^M'^ + MsM^ - M3MO + \pXM^M'^ 
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which gives the same result as that obtained through the normahzation of Isgur-Wise 
function at the zero recoil point [7] (the discussion about the extra factor appearing in 
this equation is the same as that given in the previous footnote). 
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